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ABSTRACT 



This paper is an attempt to develop a method for 
determining the slamming pressures and loads on ship hulls, 
when the relative vertical velocity relationship between 
the wave slope and keel is specified. A basic two-dimensional 
theory for hydrodynamic impact is utilized which avoids 
linearization of the free surface. The mathematical model 
for the theoretical developement is a wedge. A method of 
analysis is presented which yields the free surface shape 
and pressures associated with symmetrical penetration of a 
fluid surface by wedges of varying deadrlse angle. The 
anaylsis is based on the principles of continuity of fluid 
flow and similarity. The principle of continuity makes use 
of the fact that the volume of displaced water above t{ie 
original water line. The principle of similarity accounts 
for the proportional expansion of the free surface with con- 
tinued penetration of the body. The free surface (consisting 
of both wave and spray) is considered to be increasing in 
such a manner that the slope and velocity of the surface is 
determined by the position and penetration velocity of the 
impacting body. Once the correct free surface shape has 
been determined by a graphical iterative method, for any 
instant of time, the associated velocity potential and 
pressure can be determined numerically for various points 
on the body. 

Since only V-shaped ship sections can be approximated by 
wedges, a method is developed by which the pressures on a 
D-shaped ship section of a low deadrise angle can be calculated. 
Results are presented which resulted from calculations 
carried out on a ship with low deadrise angles in the foward 
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sections. For comparison, a calculation was performed 
for a Liberty ship section which could be satisfactorily 
approximated by a straight wedge shape. Also, the Integrated 
pressure load was determined for several Instants of 
time for the ship being considered. 

The results of the calculations seemed to give 
realistic values for the pressure distribution. However, 
the numerical calculations, as developed in this paper, 
are quite tedious. Computer solutions need to be developed. 
Also, the approximation of ship sections of low deadrise 
angles by blunt nose wedges is questionable, since the 
spray root thickness has to be approximated. Therefore, 
the free surface shape in the spray root area Is an 
approximation. The theoretical results should be closely 
checked by experiments. 
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NOTATION 
Points on the free surface 

Nominal wetted half width of the impacting body (ft. ) 

Acceleration due to gravity (ft/ sec. > 

Wave height (double amplitude) (ft. ) 

Normal draft (ft. ) 

Ship length (ft, > 

Origin of coordinate systems 

Edge of sprays origin for a 

Pressure at any point In flow field (lb/ In ) 

Distance from origin to point In fluid (ft. ) 

Wave amplitude (ft. ) 

Amplitude of complex number 

Distance along free surface curve from it® origin 
at the solid boundary (ft. ) 

Time (sec. ) 

Period of encounter between waves and ©hip 
Velocity component in x- direction (ft/ sec) 

R esultant velocity at any point in fluid 
s‘/ v 6 + w " (ft/ sec) 

Velocity component in y '•direction (Si/ sec) 

Ship speed (ft/ sec) 

Water surface velocity relative to model bottom 
(or flat plate); a function of y > c (ft/ sec) 

Velocity component in % -direction (ft/ sec) 

Linear displacements of center of gravity of ship along 
fixed axes 

Coordinate system fixed in ship 
Heave amplitude 
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Relative position of ship keel with respect to water 
surface {ft) 

Relative velocity of keel with respect to water surface 
penetration velocity {ft/ sec) 

Angle between velocity vector » U, and radius r 

Deadrise angle 

Angle between radial velocity and fee surface velocity 
past fee radius 

Spray root thickness; phase angle for heave 
Phase angle for pitch 
s -coordinate of particular fluid particle 
Wave length {ft) 

y -coordinate of & particular fluid particle 

Maas density of fluid { ) 

ft 

Pitch angle 

Pitch amplitude 

Frequency of encounter between the ship and waves 

{rad/ sec) 
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I. INTRODUCTION 



In recent years there has been a distinct trend toward increasing 
the speed and displacement ol ships. Quite often a ship's speed is not 
governed by its installed horsepower, but by its response to the sea- 
way. The term ’response ' includes not only the ship's motion, but al- 
so the effect on the ship's structure caused by the motion. One possible 
result of a ship's motion in a seaway is the slamming fore© exerted on 
the forward portion of the hull structure. The force is usually caused 
by the emergence of the bow and subsequent impact with the water on 
reentry. The ©hip may vibrate for some time after such an impact, 
and plates on the bottom shell may be damaged. If the impacts are 
severe enough, serious structural damage may occur, and even if the 
blows are of small magnitude, but frequent, the hull will possibly be 
weakened by fatigue. 

The earliest theoretical work on slamming loads resulted from 
studies on the impacts sustained by seaplanes when they landed. The 
quantitative evaluation of the impact forces was baaed on the ' expand- 
ing plate * and "spray root ' theories. The expanding plate theory origi- 
nated with Von K arm an la 1923 (1) * By this theory, the water flow 
pattern around a wedge penetrating the water surface is taken to be equi- 
valent at each instant to the flow about a flat plate of the same width as 
the wetted width of the wedge. From this analysis, the velocities, accel- 
erations, and resultant pressures and hydrodynamic reactions are ob- 
tained. Apparently Von Kerman's theory is good for low deadrise angles, 
but is only a crude approximation at higher angles <2). By the spray root 
theory, developed by Wagner in 1932, the local flow phenomena oc curing 

*) ‘ “ ' ‘ 1 * 

Numbers in brackets refer to bibliography. 
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at the edge of the wetted area is described. The high pressures that 
exist near the edge of the wetted area during impact are connected 
with th© formation of a spray Jet in which the water is quickly accel- 
erated to a high velocity. 

Seaplane impact theories have been adapted to the problems of 
•hip slamming (3, 4) in an attempt to pre-diet the pressure response 
of a given hull form. Experiments have been devised and carried out 
(4, 5, 6) utilising ship models. Th© experimentations have been carried 
out under conditions allowing three dimensional flow about the ship hull, 
which are obviously the conditions met by ships at sea. However, three 
dimensional flow analysis for a general curved body is much more 
complex than that of two dimensional flow. Therefore* theoretical 
work ha a been primarily concerned with the flow about two dimensional 
ahip sections of finite length. The slamming loads for the individual 
section© have been combined by strip theory to obtain the pressure 
distribution on th© total area of the ship bottom affected by the slam. 

There is a pressing need for the extension of the available two- 
dimensional water entry theories, in the expanding plate theory, the 
free surface is represented by a sera potential line. The impact of 
ship hull shape bodies on a water surface will usually result in the 
formation of a curved free surface, and a spray root. Theoretically 
the presence of a highly curved free surface may not b® represented 
by a straight sero potential line. It seems likely that the free surface 
is nearly a horisontal Une only for bodies of high deadrlse angles, and 
that the spray root generated at low deadrlse angles disputes the as sump 
tion of a horisontal aero potential line. Ir. general, for any of th© fit- 
ting ! techniques (ellipse and circle fitting a© well as the expanding 
plate), the agreement between theory and experiment has been rather 
poor (2), even for simple wedge shaped bodies. Therefore, a more 
exact solution of the hydrodynamic impact problem would be desirable. 
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The object of this work la to formulate a basic two-dimensional 
theory for hydrodynamic impact which avoid® linearization of the free- 
surfaee boundary condition, and which can be adapted to the ship im- 
pact problem. The theory is based largely on the work of Pierson (7) 
who developed an iterative numerical solution for a wedge entering 
the water at constant velocity. While a wedge ia used as the mathe- 
matical model in this paper, the velocity is not considered to be con- 
stant. 

The problem of determine slamming pressure® actually consists 
of two parts; first, the des cription of the ship motion ia a seaway which 
leads to slamming, and second, the description of the pressures and 
forces oa the MU caused by slamming. Ship motions are mt treated 
here, but an expression for ship motion is needed to determine, at any 
instant, the relative vertical velocity of the ship with respect to the 
water surface. Such an expression was obtained from experimental 
work for a Liberty type ship, discussed la reference (3). 

Because of the difficulty of analytically describing the flow about 
an arbitrary convex body, the expressions for the velocity potential de- 
rived ia the paper are based on the water entry of a wedge. The dead- 
rise angle of the wedge Is taken to be approximately the same as the 
dcadrUe angle of the ship section under consideration. While the above 
assumption appears to be somewhat gross, in the light of past experi- 
mental and theoretical work it does not seem to ba bad. Both experi- 
mental and theoretical calculations hare shown that the highest trace - 
verse pressure distributions occur on a giver* ship section shortly after 
it hits the water (S, 6, 8). The theoretical calculations were based oa 
the expanding plate theory, with a modification tor the hull section cur- 
vature and added mass of water. However, the portion of the hull that 
experiences the largest pressure loading is the bottom portion which. 
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for most ships* is roughly equiraimi to a wedge of low deadr-ts* aa^U, 
or can at l$mi b* approximated (at least for the period just before and 
after slamming) by a w edge of a specified angle. It is not sufficient 
just to compute the pressure distribution at the instant of slamming, 
but pressure distributions for several instant® before and after slam- 
ming are also required for a complete picture of the pressure loading. 
From a structural viewpoint, it ia the pressure loading tod its time 
variation which are important, and net peak pressures which last only 
for a few micro -seconds. 

The development of the theory is presented in Chapter II. The re- 
sults of th« theory, as applied to a given ship for which other theoretical 
end experimental data are available, are presented tod discussed in 

Chapters IE and IV. 
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H, PROCEDURE 



A, General The theoretical calculation® are carried cut for a ship 
approximately the else of a Liberty ship. The actual ship lines were 
calculated from the series of sections far which Professor F. M. Lewis 
obtained hie values of inertia coefficients <S). Professor J. £„ Kerwin 
obtained the offsets and polynomial egressions for the ship lines by 
the use of the TOSG IBM computer, Also, a ealc elation is performed 
for comparison for a section of tha Liberty ship used in references 

( 3 * 4 ). 

B. Deter min a ti on o f Ship Motion The coordinate system used for the 
determination of ship motion i& the same one used In reference <18)., 

The coordinate origin is taken to tee at midships. The <x,y,s) coordi- 
nates represent linear displacements of the center of gravity along 
fixed axes. In this paper, only heaving and pitching motions are con- 
sidered, therefore only the s coordinate is of interest. The a axis is 
vertical and t is positive upward. Of the three angular displacements 

< roll, pitch, and yaw), only the pitch angle is of interest. Positive 
indicates bow up. Also, the (X, Y, 2) coordinate system, which is re- 
lated to the ship’s principal axes, is introduced. The Z -coordinate is 
positive up, the X -coordinate positive forward, and the Y -coordinate 
positive to port. The pertinent coordinates and their positive directions 
are shown in Fig. 1. 

In order to obtain .a complete picture of the ship motion at any 
instant of lime. It is accessary to know the instantaneous draft and re- 
lative velocity between the keel and water surface for my station. 

The pitching motion is assumed to b® 

y * y cos {CO ? t + € ) 
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FIGURE I PRINCIPAL COORDINATES AND AXES OF 
REFERENCE 




FIGURE II RELATIVE ELEVATION OF KEEL ABOVE 
DISTURBED WATER SURFACE 
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and the heave is assumed to be 

* « E COS {CJ a t + S ). 
m 

The wave elevation at any station along the ship is given by 

T T OOS<~p- rUJei) 

where oj^ * * (Y+ ) » period of encounter between ship and waves 

in head seas. 

The relative position of the keel above the disturbed water surface 
is given by: 

% *r- S - XV + H (1) 

x 

for may station (see Figure II). Emergence is associated with a < 0. 

The relative velocity between the ship and wave at & distance X 
from amidships is given by 

* * T i-X? (2) 

where T* * ~~ ° - u) 0 *f sin (u>ci + ), 

o w m A 

and Cd 0 is the absolute angular frequency of the wave and equals , 
and * « • - * r sin ;a; e t -5- » ), 

and [y * ~jr " - 4^ cc>e sin (&V +6). 

From reference ( 1 1), experimental values of fi and were obtained 
for a 5. & ft. model of a Liberty ship in regular waves of a specified 
height. The values given for the model are the same for the full scale 
ship, with the exception of the heave amplitude, which is scaled up for 
the ship by a ratio of the lengths. The pertinent variables for the mo- 
del, full scale ship, and equivalent Lewis ship ore tabulated in Table I. 

The calculations of the relative velocity between the keel and the 
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water by £q. (2) will result in a sinusoidal variation of velocity for 
any station. While this is not correct* it is shown in reference <3) 
that this is a good approximation for the relative velocity of a ship 
subject to slamming motion in regular waves. In other words* the 
change In relative velocity due to slamming loads is very small. 

Slamming, as generally defined, is the sudden change of the 
acceleration of the ship. While the relative velocity does not deviate 
much from that predicted by (2), the influence of slamming on the accel- 
eration is quite marked (3). However, in this study of slamming, the 
pressures are governed by the velocity and penetration of the keel in- 
stead of due acceleration. For the ship hulls under consideration, the 
portion'', of the ship between stations 0 and 5 is subject to slamming. 

From reference (3) the instant of slamming ia determined to be in the 
vicinity of t * 1. 20 sec. The penetration, and relative velocities of sta- 
tions 0 through 5 are graphed for the interval from t » 1.20 sec. to 
t » 1.44 sec. in Figs. HI and IV. These values will be used to deter- 
mine the slamming pressures on the forward portions of the ship. The 
values used to plot the curves in Figures III and IV were computed from 
equations (la) and (lb) which resulted when the amplitudes and phase 
angles from Table I were substituted into equations (1) and (2). 

2*x 

% x * -12.8 cos {co e t + ) - 8. 38 cos -0. 623)-0. 122Xcos{^t+0. 25 

(la) 

2-v 

* * -8.3 sin(cygJt+ ^2~ sin (<Ve<t - . 823)- . 122 Xsln (cv e i + . 252) 

(lb) 



C. Theoretical Analysis 

Conditions of the Problem 

1. The problem is treated as a two-dimensional one. The velocity po- 
tential and pressure distribution will be calculated for two-dimensional 
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TABLE I 





Model 


Liberty 
. . ship 


Lewis 
type ship 


Length (BP) L, ft. 


5.5 


416.0 


420. 0 


Beam, B, ft. 


0. 732 


56.9 


57. 5 


Draft, H, ft. 


0. 21 


15.9 


16. 0 


Scale factor 


1. 0 


75.7 


76. 5 


Speed V, knots 


1. 1 


9. 6 


i 

9.7 


Wave length X, ft. 


5.0 


379. 0 


382.0 


Wave height h, ft. 


0.334 


25.3 


25. 5 


Pitch amplitude , degrees 

1 m 


7 


7 


7 


Heave amplitude Z^, ft. 


0. 007 


5. 30 


5. 35 


Pitch phase , degrees 


14. 44 


14.44 


14.44 


Heave phase & , degrees 


-35. 7 


-35. 7 


-35.7 


Frequency of encounter 6U?, 


8. 74 


1. 01 


1. 00 


Period of encounter T , sec. 

e 


0. 72 


6. 28 


6. 28 


rad. 

Wave frequency , gec 


6. 36 


0. 73 


0. 73 
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sections (each a station spacing la length), and then combined by strip 
theory to yield the pressure load on the forefoot, of tho ship. 

2) Each impacting section is assumed to b® perfectly rigid and to be 
symmetrical about the vertical centerline plane. 

3) Th® velocity of penetration is not constant, bat varies according 
to equation (2). 

4) The flow field is assumed to be irrotational (a necessary condition 
for the existence of a velocity potential). 

5) Th® fluid is assumed ia be frictioaless and incompressible, and 
the velocity of penetration is assumed to be high enough for the effects 
of gravity and surface tension to be neglected. 

6) On the basis of the assumption of irrotatioaal flow, the particles 
of the original surface remain on the surface and no new particles are 
added to the surface (12). Since the fluid ia incompressible, the arc 
length along the surface from, the impacting body to any given particle 
is constant and independent of time. 

The preceding conditions do not define the shape of the free sur- 
face. However, for a wedge shaped body, which ia the mathematical 
model used for the theory, two general conditions can be formulated 
which will define th® shape of to* fre^ surface; 

a) Continuity : The fluid displaced by the body must appear above 

the original water surface in the form of a wave and/or spray. The arc 
length of the surface must be constant. 

b) Similarity; Since the impacting body is symmetrical, and th® 
immersion is normal to the surface, the field of flow is divided into two 
halves by the centerline of symmetry. By neglecting the gravity forces, 
the surface shape and dynamic state of the field of flow can b® defined 
entirely in terms of the penetration and velocity of the body (?). If the 
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penetration velocity ia constant, the flow field and surface shape must 
be geometrically and dynamically similar at all time®. However, in 
the present case, the penetration velocity is not constant, and there- 
fore the surface shape is assumed to be only geometrically similar at 
different time®. 



Application of General Conditions 

CONTINUITY The general shape of the free surface due to penetra- 
tion by a wedge shaped body is shown in Fig. V. A* shown in the figure, 
the free surface can be thought of as being made up of spray and a wave. 
The two merge together in the ’spray root*' area. Ia determining the free 
surface shape, the wave rise i$) is given by equations 11 and IS of refer- 
ence (13). 



2 c tan (5 






(3) 



This is the same expression which gives the surface rise in the ' expand- 
ing plate" theory. The wetted hall width (c> is the point on the body which 

is wetted by the piled up water in the absence of spray. The wetted half 

% c 

width at the original water level is ----- (Wagner* see reference 14). 
The continuity condition is satisfied if Area A * Area B, and if a * y. 
However, the continuity condition Itself is not sufficient to determine 
the free surface shape since any thickness? of spray may be associated 
with an unlimited variety of waves, and still not violate the continuity 
requirement. Also, the free surface in the spray root area, where the 
surface wave and spray merge together, could assume any one of many 
shapes. 



Similarity The principle of similarity for the proportional increase of 
the wave and spray with continued penetration is the second considera- 
tion. Not only must the displaced fluid appear above the original water 
surface, but the wave and spray must be increasing in such a maimer 
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FIGURE V FREE SURFACE DUE TO PENETRATION 
OF A WEDGE 
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that the shape sad velocities are everywhere determined by the position 
sod penetration velocity of the entering body. By a geometrical con- 
struction and integration process proposed by Wagner (see reference 7), 
it is possible to relate the surface shape and velocities,, so that for any 
spray thickness, one corresponding wave and spray-root shape can be 
determined which meets the requirements of continuity and similarity. 
However, Wagner's proposal was for a constant entry velocity, the follow- 
ing theory is developed for a variable penetration velocity. 

In m ideal fluid where gravity, viscosity, and surface tension are 
neglected, the shape of the free surface is uniquely determined by the 

shape of the impacting foody, and the siae of the shapes is determined 

a 

by the amount of penetration . 

Consider the several stages of the penetration process of a wedge 
shaped body Figure VI, using an y-s coordinate system. The penetra- 
tion of the body into the fluid causes the initially flat surface to rise as 
a wave, and the nearest particles are deflected out along the wedge sur- 
face as spray. The resulting surface shape expands as the wedge pene- 
trates deeper. The resulting surface shapes are o\ ajfoi, o’* a? b*, etc. 
For the same relative position In the fluid and on the surface, the pres- 
sure and velocity depend upas. the surface shape, and mast bear the same 
relationship to the penetration velocity at any instant of time (1), The 
above statement applies to the relative position in the fluid, and not to 
a particular fluid particle. The positions are indicated by o'* , o'#, o'® 
and as , a* , a*, etc. If the penetration velocity was constant for the body, 
the velocity of the fluid at corresponding points would be equal. That is 

e) ’ " " " ’ ‘ “ 

Of course this statement is only true for a certain definite time inter- 
val after impact. Ae the body penetrates deeper, much of the free sur- 
face shape will be dissipated a® spray, and as the penetration velocity 
decreases, the influence) of gravity will be felt. 
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U * U * U„ * K * where K represents a proportionality factor 

between the velocity U at positions a t , a 3 , etc., and the penetration 

velocity, s . However, for varying velocity, U j U f U , but the 
x aj a? a 3 

velocity at corresponding points will still be proportioaai to the pene- 
trating velocity. 



However, it is not unreasonable to assume that K * K a K , K. * 

sli a * a / 1 >| 

K.» K. , etc. Therefore, the velocity at corresponding points will be 



directly proportional to the instantaneous penetrating velocity. 

The velocity at any point on the free surface, such as a® or b s , 
is determined from consideration of this condition of similarity. Con- 
sider point to* for example, A particle on the surface at point b 8 is 
instantaneously moving along the radius with a velocity and at the 
same time sliding past the radius toward the spray region with a velocity 
U gfe (see Fig. VU). 

The vector sum of the two motions <U + U ) yields the instan- 



taneous resultant velocity of the point in question. 

The radial expansion velocity, U , the velocity of the surface 

T 

past the radius, U , and the resultant velocity, U, are each directly 
proportional to the penetration velocity, z . Therefore, for constant 
entry velocity, or for an average velocity taken over a small time in- 
terval, the velocity component® (U , U ) at a position are constant, and 

T $ 

the vector sum may be obtained from the radial distance and the arc 

length. 

Consider the point b*, on the curve in Figure VI. The radial 
distance r is given by: 




u b,* K b i x <t * ) 



( 4 ) 



**fe «b' 
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The vector s is constructed tangent to the surface at b a so that 
s equals the arc length 0'b ; ^ 



H 




dt 



The vector sum is then: 

a +■ r * 



J 0 



( 5 ) 



The individual fluid particles do not follow the radiating lines, 
but remain at a fixed are length from 0* . The motion ef a particle is 
therefore a function of its relative position on the free surface, and in 
given by the following Integral relation® 

§ -&> */v4i (0a) 



V 




<6b) 



Equations (0) are referred to the y and z axes in Figure VI. The motion 
in the y direction In a time interval t is expressed by f -£>,£> being 
the original position of the particle. The symbol q represents the vert- 
ical elevation of a particle (a -direction). Mow the velocity of a particle 
(v, w) Is to be determined for any point in its path, so that Equations 
(6) can be solved. There does not seem to b© a simple analytic solution, 
but a graphical solution can easily be carried out. However, the solution 
can be simplified by changing the variable of integration from time to arc 
length. This is done in the following manner. If the penetration velocity 
was constant, the arc length and times at two corresponding points on 
the free surface could be related by: 



S, 






°sb“ 



k 

ti 



s 

h 



(reference ( 7 ) ). However, for varying velocity at time t, H ^ * K g ^ a x Ui ), 
where is the proportionality factor which relates the velocity of the 



- 10 



S A 



%f v Vrr% i »* ». - v^.^. SiTT 

• ' -,^>i -f ■>.»* *iao/y,» % 

*** 



. ■*> * ' ■ • fi •»» ♦<fT 



tu •• ji • t • i 

Cr- 

,M^ i'. ■ • • * '.V. Ji . - • 1 I «.« >4- U. Mf lV*J VJT 

•j ?i 3*>.*45 * V-. : I - , •..,*.>» f - • - i -d 

• -• * ' 1 ' ••„• y, 'r .. . ,■! • 



Uii 



*.d £• 



I J > U ifAli . d* , i u*Viy 



.V l 



.. . *■ : J\ • r ; ...» ■ . » •• . " ■ * <*) 

»•’ "» ; <S. cl ^ ^ W* >np< * I tow •: i »«II * jJ n V. ' >/-l Cl 



•rw %,* **»< * lo4i.;» *<i X , air*!**) «ci ^ iuu^iiq *di 

1 1 -■ . ■ •■ .*■.-.»*. stto wv* '■*. J - . ••-*.•: *» • • * >r. 4 i»vf 4 t toal 

-M'iJh . * ■ r»*Q •< . . ,* ; f ‘w ,*r) 



1 •■ a » *j i , , j' : . , < .. • ju ojl (S) 

0~AiLte i .»*■ < 5 -* <~m # • - Mi #• • - I ... A * t*/d 

«-» X> tflC a J*- »•> ** * ■ t 5 , ■■• #d <L «0 

\ 1 (?c 4 »v *(« V - y '* il' •'• :» ♦ *#&■ i . it . . di^r.ry | 

it* nn;' ,« i*ta? >rf» it %0t in o|.-»i ra 4’ . >«a».al9 j*w 

t .**1 1 •STkl *fl 5 



•I 4 ^i ^ <9 Mrir to fltootoi i mm v 4 t wwh > .1 iv> ttwitin) 

W 0 ll |MMIi «to Mitoi <*• 10 MI itUlNltoMfnq Ml to . | MMftl 



fiurfafs sliding past the radius Ob to penetration velocity. However, in 
the time interval from U to t s , we can assume that U . is constant and 

30 

equal to its average value. According to reference (IS), aat average valve 
for a function, f{£), can be defined sm follows: 



y° f (t) * 
t> + U 



f (t+ ~ ) + f (t - -~) 



Therefore, at time —~z — , th« average value of U . 

& SD 

to be constant over the interval i A to . For any time. 



% m u <*L ± ii. ) * ?sl 
t 3 sb 1 2 1 t 



can be assumed 

ta , 



By this reasoiiing, an expression for arc length in terms of time can 
be obtained for any time interval., At : 



where t * the earlier time, and S * corresponding arc length on the snr 

face at time t a . 

St, 

Also, dt. * - d s 



Substitution for dt in Equations (6a) and (6b) yields 




Equations (fa) and (7b) can be solved numerically m follows: 



(7a) 

(7b) 



Referring to Fig. VI, it is seen that the paaetr&tion velocity has 
approximately a linear variation over the time interval under considera- 
tion. Since the particle velocities at corresponding points are directly 
proportional to the penetration velocity, they can also be assumed to 
vary linearly. Using this assumption, the integral in Equation (5) can 
be expressed as 
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K dt * V 

where is the average value of U, o for the time interval under question. 
From this, the average value of the velocity components, in Equations 
(7a) and (7b) cam be found (vt and wt). Therefore, a step by step inte- 
gration of Equations (7a) and (7b) can be carried out as shown in Appen- 
dix A. Any trial free surface shape which satisfies continuity can be 
used for the determination of the surface velocities. The velocities 
are then used to determine a new surface shape. This procedure is 
carried out until there is a negligible difference between assumed and 
computed values of surface shape. 

The integral equations (7a) and (7b) relate the surface rise at a 
point to all the points more distant from the body. Therefore, the path 
of each particle on the surface follows the same related path, and the 
conditions for similarity are satisfied. It would appear that the integral 
equations reduce the number of spray-wave combinations to one for each 
spray thickness (7). 

The spray root thickness itself was shown in reference (7) to be 
related to wetted width (c) and deadrise angle (j3): 

6 x c ( ~ )‘ (8) 



Potential on Free Surface 

Referring to Figure VII, the tangential velocity, along the 
surface at a position, b is given by - ® _ + IJ^ cos Y . For points along 
the surface (at a specific depth) cos Y = ~ * . 



Therefore, 



d S 



a 



U . - U . cos Y 
sb rb 



U , - U 
sb 



rb 



dr 

ds 



(9) 



*) dr 



a 



ds 



for a particular surface should not be confused with 
constant position, such as b. 



s 



r /s for 
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Since and °rb ere proportional to c^ a and k^ can be assumed equal 
to a mean value over the time interval under consideration, U g and U 
can be considered approximately constant over the same time interval. 
A differential equation can be obtained by multiplying Equation (9) by 

tda. UUs^s, t U cerh 

T 

- td« « sds - rdr 

Integration yields 



which is an approximate expression for the velocity potential on the free 
surface with varying penetration velocity. 

Potential and Velocity Distributions on Body Surface 

Assuming that the free surface shape for a given deadrlse angle 
is unique, and is a result of potential flow, the velocity and potential at 
any point can be represented by a complex analytic functions 



The Cauchy integral equation may fee applied to the complex func- 
tion m described in refox s ence (?), By Cauchy’s integral formula, a value 
of an analytic function in a region is determined by its values on the boun- 
dary. There is no choice of ways in which the function can be defined 
at points away from the boundary once the function is defined on the boun- 
dary. Every change of values of the function ai interior points must be ac- 
companied by a change of its values on the boundary if the function is to 
remain analytic <18>„ The integral formula for the point within the boun- 
dary la; 
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m is the complex coordinate oi the position in the field, Sine* point t- on 
til* solid boundary are being investigated, lb® bounding caniour circles 
only half way around the point so that iSsr* becomes ari, Also, f(tt) °> u t- iv 
* IK*) on the boundary, and equation < 12) becomes; 



0 <* > - 

o 




(IS) 



i c 

The expression can be separated into the real and imaginary part® by 

placing the origin of eoor&nafrfta at the point & under consideration, 

o 

and writing * in the polar form; 

i& , in A + ie 



a * He" * «' 

m*J* m * ie 



[d O oE) ♦ id*?] 



(14) 

US) 



Substitution of the expression for * and da into equation ( 5.3) yields 

U (s ) • V - iv * f u - ir fd (lnA) ♦ id&l (16) 
q o o al L J 

Rearranging tarsaa give* 

itru 4 ffr * J"d (laB) 4 vdS 4 i J^ud^ - vd UnE)J , 
so that tie real and imaginary ports may be separated; 

v ; * ~ J vdd 4 Ud (InB) (l?a) 

u * J' d£ -vd(lnR) (l 7b) 

The construction for performing the numerical integration of the 
above equations is shown in Figure VXTX The complex coordinates were 
oriented so that the y-aicia (imaginary axis) is al^rg the w©d&e, positive 
toward the spray tip. The x~n*it (real axis) ir. norm®! to the body and 
positive oat, 'The potential on the body surface may fe* obtained by inte- 
grating the velocity along the wedge from the apray tip to the apex. 



Pressure Distribution 



The general equation for fluid pressure in a potential flow of in- 
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FIGURE VIII POLAR COORDINATE SYSTEM FOR THE 
DETERMINATION OF THE VELOCITY 
POTENTIAL ALONG THE BODY SURFACE 
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compressible fluid, in used to d«termiu< the pressure distribution 
oa the foody ( 1%) 



9 





+ F it) 



(18) 



The gravity force <c ) wd extraneous impulse# (F (t) ) are aero since 

gravity has been neglected, and the fluid was assumed to be at rest at 

infinity. Therefor* Equation ( 18) become# 

If 
2 * 



P* !i 
f t 



U9) 



The resultant field velocity (U) can be determined on the free surface 
(Equation (8) ), and on the foody (Equation (17) h Now it remains to de- 
termine an expression for —• . 

-* v 



Assume that the depth of penetration at time t x is given by the aver- 
age penetration velocity (* ), in the interval t * o to t ■* t 4 , multiplied 

% 

by t* . This is the same as assuming constant velocity penetration in 

fgtj 

the specified time interval. At any position, {y£ f in the fluid, the 
potential is therefor© a linear function of time in the interval to i% (7). 



For any time interval t* to t, assuming s^, 

4*1 * 17 <20) 

However, for a fixed point <y„ s) in the flow field, the above expression 

is not correct. If the point <y, s) coinrides with the position r. , * at 

* i 

tj * at time t ■ a t the position r^ + *^j 1® a dUtaacM Ar * — ^ — 

from (y f #}. For ia»U A t, the resultant vector velocity caa foe asst ned 

constant in the region. 



At t * ix $<y, a) •* and at t * t% + a* 4^y, ~ ^ CORa ^ r ^ 



*) 

See Figure XX, <T\ a ) refers to the distance to the same relative point 
In the fluid at a deptli* oi° penetration a , Therefore, as time increases, Tl 
increases. 
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FIGURE IX FIXED AND RELATIVE POSITION POINTS 
IN THE FLUID 
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U cos a is the velocity component along . Therefore, 



Mf» ») t, ( f 4t *> (4 • 4^ - cosa *• t (21) 

From Equation (20) 



s* ta 

Equation (21) then becomes 

f« .. 



* ti 4- At x At . 

Tt * ~ A “ I — ** 4^, + ~r* *. 



^ U ft, 



At * ' r •* Uhi co» a 

n 7 t» 

Taking the limit of Equation (22) as At approaches Or 

si . liM &L. . "* < <ft . 

At-»0 At i t-40 r — 

ti 



( 22 ) 



&r 



0 cos a ) 



£«■ a t 

>t v T 
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U cor- a 



(S3) 



Equation (S3) now applies to all point® in the fluid. The pressure distribution 
is given by substituting Equation (S3) Isa Equation (10). 

r u s 

“ - “jr + ~£~ U cos a - -g- (24) 

On the body surface, the pressure may be found by using the values of , 

4, U, and a already determined. The total load on the body is obtained 
from the integration of the pressure across the projected wetted width. 



Applicati o n to Sh ip Sect ions 

There are two possible methods of solving the water entry problem 
of bodies ©f general shape by considering the free surface shape. The 
first is to construct the free surface based on the boundary of the actual 
ship section as it enters the water. However, this method is questionable 
since the similarity principle has not beers proved to hold for a body of 
general shape. The second method is to approximate the ship section by 
a wedge-ilk® shape of m equivalent dc&dria© angle. 

The period ©f Interest In a a lamming investigation, for a certain 
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ship section, ia the time period immediately after impact with the 
water. The maximum impact fore® is produced either at the instant 
of impact, or immediately after entry. For the Liberty ship investi- 
gated in reference (3), it was found that the slamming pressures were 
of little significance by the time the draft at station 6 wm 9, S ft. Ia 
fact, the time Interval of significant pressures was less than 0. 17 sec. 
For this problem, the time interval under investigation is 0, 24 sec. 

For a hull section with sharp curvature, the spray will not neces- 
sarily adhere to the boundary. From photographs of wedge impacts ia 
references <14, 17), the spray, after it leaves the wedge surface, shoots 
off at an angle roughly equivalent to the deadrise angle (see Fig. X), 
Therefore, it seems possible that V-shaped sections, or sections with 
fine lines, could be approximated by a wedge of a certain deadriee angle. 
However, slamming is more dangerous for sections with email deadrise 
angle®. If Is obvious that these sections could not satisfactorily be re- 
presented by a V-shaped wedge. 

For U-shaped sections with small dead-rise angle©, the method 
for constructing the free surface *shape will have to b© modified. It is 
obvious that the wave rise cannot be estimated correctly fey Equation (3), 
Figure XI is a sketch of a typical U-shaped section. The pressure distri- 
bution for such a section is satisfactorily predicted by the expanding plate 
theory, except in the spray root region, from y * to to y * c. It is there- 
fore reasonable to expect that the pressure on the flat portion of the hull 
can be given by assuming the pressure distribution for a flat plate from 
v * o to y * b. The flow past this "flat plate” can be used to estimate the 
wave height at the free surface. From (18) the distribution of the normal 
velocity for a flat plate of width 2a, at points y > a, is given by 
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FIGURE X IMPACTING WEDGE 
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FIGURE XI FLOW NOMENCLATURE 
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V will be the velocity of the free surface particles, The were height 
At A point y can be found bys 




For the wave height determiaatl? as, the dimension a i® taken to be 
equal to c% the width of the section immersed at the original water 
line (see Fig. XI). When constructing the free surface, the spray is 
assumed to s2ioot off at an angle determined by the slope of the hull in 
the spray root area. This is the same as assuming that the side of the 
ship is straight, and that the impacting body is roughly equivalent to a 
blunt nosed wedge. Assuming that the similarity principle still holds, 
the free surface is constructed by using the arc length and continuity 
conditions, and checking by Equation (0). 

For a U-shaped section, ike velocity potential for the flat portion 
of th# hall (from y * 0 to y * b) cm be assumed to be that for irrotatioasl 
flow about a flat plate. From (it) the complex velocity potential past a 
plate perpendicular to the flow is given by; 

<*■ « - i v \J jr - a* " - \j a y - y® (Sit) 

Equation (1?) can be used to compute th# velocity potential on the hull 
from y * b to y » e. The two velocity potentials should agree at y ° b 
for any specific time. (A sample calculation is discussed la Appendix A. ) 

The calculations, from which the results la this paper were de- 
termined, were carried out for the Lewis type ship discussed previously 
(set Figure XU). 
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FIGURE XII LEWIS SECTION SHAPES 
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Included In this section ere the curves of the theoretical pressure 
distribution for the ship under consideration. These curves arc shown 
In H gures XIII * XV . Also, one transverse section of the Liberty ship, 
used as a model in reference <#), was investigated. The pressure distri- 
bution for this section is shown in Figure XVII. This section was quits 
V -shaped and had a deadrls# angle of & 80* . Straight wsdg® theory was 

used for calculation of the pressure distribution, end the results are com- 
pared with the results given in reference (3) for the same section. 

A s amply calculation for station 3 of the Lewis section ship is dis- 
cussed in Appendix A. The tangential velocity (v^) on the flat portion of 
the hull was calculated by two methods (Table IV). The first, method was 
by the assumption of flow over a flat plate. This is the method used for 
all of the low d eadrise sections under consideration. However, for com- 
parison, the tangential velocity on the bottom of the hull was calculated 
by use of the Cauchy integral formula and equation (1?). The results 
check closely with those from the flat, plat® flow assumption. 

The pressures load on the forefoot was determined by strip theory 
for several instants of time. The integrated load distribution is shown 

in Figure XVI. 
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IV DISCUSSION OF RESULTS 



A. Introduction 



Before discussing the results, the theoretical assumptions made 

will be restated. 

1) Only two dimensional sections were considered. 

2) Each impacting section was assumed to be perfectly rigid, and to 
be symmetrical about the centerline plane. 

3) While the velocity of penetration was not constant, for the small time 
interval under consideration, it was considered to vary linearly and aver 
age values were used. 

4) The flow field is assumed to be irrotational. 

5) The fluid is considered to be frictionless and incompressible, and 
the effects of gravity and surface tension are neglected. 

6) On the basis of irrotational flow, the particles of the origins! surface 
are considered to remain in the free surface. The surface deflected by 
the impacting body is assumed to have a constant arc length. 

7) The theoretical procedure was based on a mathematical model of a 
wedge with constant entry velocity. The theory was applied to the water 
entry of actual ship sections by assuming that the impacting sections 
were approximately blunt nose wedges. That is, the free surface formed 
during the initial impact stages was considered to b© equivalent to that 
formed by a wedge of an equivalent deadrise angle. The spray was con- 
sidered to shoot off at an angle determined by the slope of the hull at the 
wetted width. For the hull sections where the deadrise angle was approxi 
mately aero, the velocity potential was determined by the expression for 
flow along a flat plate. 
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8) The ship motion expressions were determined experimentally for a 
Liberty ship of approximately the same dimensions as the ship found 
from Lewis sections. 

B. Determination of Pressure Distribution s 

The ship motion expression for the relative velocity did not include 
the effect of ship forward speed. For a small pitch angle, the total re- 
lative normal velocity of the keel is; 

^ Uf + xi+2-r 

where U is the effect of the forward speed on the normal velocity. 

If the pitch angle is small at the instant of slamming, U can be safely 
neglected. At t * 1.284 sec., the first time for calculations, the pitch 
angle was approximately 0.2* . Therefore, the contribution of U to the 
relative velocity was not important. 

Figures 2211 to X V represent the pressure distribution on the fore- 
foot of the ship at approximately the time of slamming. The pressures 
were computed for all stations from 1 to 5 which are immersed at the 
particular time. The pressure distributions are not really sufficient 
to give a definite value to the maximum pressure to which the ship is 
subjected. 

The pressure distribution curves are similar to the curves ob- 
tained In reference (3) for the liberty ship. While the peaks on the 
curves obtained in this paper were lower, this is probably not significant. 
The peak pressures are of extremely short duration, and are relatively 
unimportant. For example, at t * 1.270 station 3 was out of the water, 
but .0184 sec. later it had a peak pressure of 550 psi at y B 5. 8 f , and 
.010 sec. later there was a pressure of approximately 120 psi at y * 5. S’ . 
However, it does seem that the average pressure distribution at a station 
varies directly as the peal: pressure. 
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An attempt was made to determine the pressure at station 3 at 
a slightly earlier time, but the construction o? the free surface was 
extremely difficult to perform, and the author had no reason to believe 
the results were correct. It is quite possible, however, that the peak 
pressure might be even higher at station 3 at an earlier instant of time. 

Referring to Figure XVII the pressures on station % o? the Liberty 
ship were uniformly lower than those determined in reference (3) by 
the expanding plate theory. This result is probably the more accurate, 
since the ship section could accurately be represented by a wedge during 
the time of initial entry into the water. The wedge assumption would not 
have been valid at a draft of 6 ft. , or greater, but by this time the im- 
pact pressures have been greatly reduced. From the results given in 
reference (7), it is to be expected that the computation of pressure 
distributions on wedge-shaped bodies, by considering the free surface, 
will give somewhat lower pressures than those calculated by the ex- 
panding plate theory. 

The integrated a lamming load distribution in tons/ ft. of length, 
is shown in Figure XVI. At best, these curves are crude approximations, 
since each curve was determined by only S or 4 points. A much better 
method would be to compute the pressure distribution at ~ and ~ sta- 
tions in the area of slamming. This was not done in this paper since 
the offsets for fractional stations were not available. However, even 
with the points available, a good idea can be obtained as to the total 
slamming load on the forefoot of the ship at the indicated times. 
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V. CONCLUSIONS 



1) The method of analysis given in this paper is at least a reasonable 
approximation since the curves of pressure distribution are of the cor- 
rect shape. 

2) The procedure for determining the pressure distribution by con- 
sidering the free surface shape is a long and tedious process. There 
are inherent sources of error in the construction of the free surface, 
and in the graphical and numerical solutions for the pressures, even 
for a simple wedge. The procedure, when adapted to a ship form, is 
even more questionable. However, there does seem to be some merit 
in the procedure, and it is concluded that it can be refined considerably. 
The biggest improvement would be to develop a computer solution for 
the procedure, 

3) The wave rise is probably accurately given for an impacting body 
by either equation (3) or (27), depending on the shape of the body. How- 
ever, it is evident that equation (8) cannot be used to determine the spray 
root thickness for a ship section. When the free surface was constructed 
for the blunt nosed sections, the spray root thickness had to be estimated. 
This was possibly a large source of error. 

4) The determination of the free surface shape in the spray root area 
is quite critical. This is where the particle velocities are the highest, 
and even a email error in the shape results in a large error in directions 
and magnitude of particle velocities. However, even with the possible 
sources of error in the spray root region, the pressures actually obtained 
for this region were as realistic as those determined by expanding plate 
theory. 
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5) The procedure described in this paper is only valid during the 
initial impact, and a short time afterward. After the body has pene- 
trated to a sizeable depth, the effect of gravity, and spray dissipation, 
will render consideration of the free surface useless. However, the 
highest pressures will, of course, occur at the time of impact or imme- 
diately afterward, when the added mass of the body is changing most 
rapidly. 
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VI. P.ECGMM ENBATIGNS 



1) Experimental and additional theoretical work should be performed 
to determine if similarity principles can be applied to ship sections of 
general shape. Particular attention should be given to determining how 
tiie expansion of the free surface is determined by the shape of the im- 
pacting body. As a start, experiments could be devised which utilized 
impacting wedges of various deadrise angles and a high speed camera. 
The spray root thickness and the angle at which the spray leaves the 
body should be thoroughly investigated. 

2) Computer procedures should be developed for the solution of 
hydrodynamic impact problems involving a non -aero velocity potential 
free surface. 

3) Accurate theoretical and experimental methods should be developed 
which will give the total slam load acting on the forefoot of the ship at 
any Instant of time. 

4) The effects on slamming loads of hull flexibility and water com - 
presaibiilty should be investigated both theoretically and experimentally. 
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VII. APPENDIX 



A. Sample Calculations 

B, Literature C> tatlons 
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APPENDIX A 



Sample Calculations U-Shaped Section {Lewia Ship, Station 3), 

t a 1, 284 sec.) 



The first estimate of the surface shape is computed by using 
equations (27), and the continuity principle. From the plotting of the 
wave rise, the wetted width (c) is determined to be 5.95 ft. (See Figure 
A-I), The slope of the hull section at y * c is determined to be approxi- 
mately 6S® . This is the angle at which the spray is assumed to shoot 
off. The spray tip. O’ , ie located by measuring along the wave and 
wedge from y * ID 1 a distance s * y * 10 ! . (At the distance y * 15’ , 
the wave slope is quite small, so that s * y is a close approximation. ) 

A spray thickness of £ /c » .12 waa assumed. Once the spray and wave 
have been merged together, the spray tip is located more exactly by 
following the curved surface (O' is located or. the line drawn tangent 
to the hull at y « c, at an angle of ft ® 83®). The area of the displaced 
water above the original water line ia made to equal the area of the 
submerged body below the original water line. Then equations (6) are 
used to check the position of the water particles. The steps of the num- 
erical integration for the surface shape are shown in Table II, The 
new particle positions were determined by the following approximate 
forms of equations (6): 
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After the free surface shape ban been determined, the velocity 
potential can be determined on the body. For the flat portion of the 
hull, from y * 0 to y * 4. 7, the potential is determined by equation (28). 
In this computation, a is taken to be equal to the wetted width at the 
time under consideration, but the equation is only used for the flat 
portion of the hull (from y « 0 to y * 4. 7). 



For the portion of the hull (from y * 4. 7 to y * 6. 05) the velocity 
potential is computed by integrating the velocity along the hull surface. 
The tangential velocity (v q ) along the hail is determined from equation 
(17). Four points are chosen on the curved portion of the hull for the 
determination of the velocity (see Figure A-l). The iy and x axis are 
constructed tangent and perpendicular to the hall at each point under 
consideration. The velocity (v^) is determined by equation (17) by num- 
erical integration around the free surface; 



v o * vd<9 + (In??) 

The steps in the integration are shown in Table HI for point 2 (*, v, end u 
refer to the position and velocities of points on the free surface),. On the 
body, v is equal to s cos |3 at any point. Knowing the component velo- 
cities, the resultant velocity (U) can be found at any point. The velo- 
cities and velocity potential determined at the points on the hull, shown 
in Figure A -I, are listed in Table IV. 



Now the pressure can be determined for equation (24); 





r 

+ ” U cos a 







where is the distance from the origin to the point on the body under 
consideration, and a is the angle between and U. The results of the 
solution of equation (24) for station 3 at t * 1.284 are listed in Table IV. 
The time, t, involved in equation (24) is the elapsed time from impact 
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of the body. For the case under consideration, t (impact time) * i, 270, 
and the computation* were carried out from time t* * 1. 284 sec. There 
fore, t * tj - t 



* .014 sec. in equation (24). 
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TABLE II 



NUMERICAL INTEGRATION FOR SURFACE SHAPE 
Station 3 t = 1.284 z y = 0.3 z* = 26.0 @ = 6j>° 
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wt A s 
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Ut 
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8 


.17 


.60 
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.25 


85 


1 .00 


6.5 
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.018 
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.35 


79 


1 .90 


6.25 


.45 


4.00 


.103 


.026 


.080 


.50 


69 


4.20 


6.05 


.59 


5.80 


.154 


.026 


.106 


.64 


55 


6.70 


5.9 


.71 


6.50 . 


.187 


.037 


.132 


.77 


45 


8.30 


5.7 


.89 


6.00 


.184 


.038 


.169 


.96 


27.5 


10.55 


5.5 


1 .09 


6.00 


.198 




.207 


1.15 


25.3 


10.55 



Area of submerged body = 1.54 ft 



Area of displaced water = 1.60 ft 
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.53 


8.70 


.268 


.056 


.081 


.46 
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3.70 


.288 




.137 


.75 
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TABLE III 



Numerical Integration of Equation (17) 
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ud(ln R) 


s = 10 
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- 5 


.4o 




4.2 


~1 .43 












0 
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s = 8 
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- 20 


.40 
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.79 












- 8. 






33.0 


s = 6.5 
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- 42 
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- .29 












163 
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475 


25 


.84 
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1.27 
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27 






230 


s = 0 
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354 


1.57 








1 



600 910 



(point 2 ) = ^ 7 ( 600 ) + 77 ( 910 ) = 481 ft/sec 
u^> = s x cos 63 °= 11.8 ft/sec 



TABLE IV 



Velocity, Velocity Potential, and Pressure Distribution 
on Station 3 of Lewis Ship(t = 1.284 sec) 
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1 1 .8 


200 



' 50 - 








t * — 1--- 






■1 

_ . ■ 





. i 


f 
















i 


















































































































i 


i 


















i 








i 











































M i i 


i 


1 


















“7 


u 


/ 














































i 
























1 






‘ i 


- 
















































U 


!□ 








i 














I 














. 












M_ 


! 


















Jj 






I 


i 




'T 

! 














/ fpi 






1 








































MM 










imi 




















j i 






i_ 










































































i 


































| 
























, 




























I 








_ 




























j 




































I 
























































































M 


















_ 






















; 










| 
























~r 














M 




















„ 




































































LD 








i 














































































" . 


















































1 




















"i 










-J 




















































d: 






















1 , 






r 




' 1 : 








1 
























M 






























1 




[Ml! 


ZT 1 

1 u 








1 








































1 


J 


— \ 




. 






! | 
































! 


i 




ID 








j i i 


n 








1 




J 


£ 


■ K 
O 


































1 


[ 








; 1 


i 




1 








\ 


i 1 i | 








f | 
















: 
























' i 






<j£d 








i 




























i i ■ 










| ; f i | 








X 




of 

M 


-7 - 




I 






M M 1 


i 














i-i -i i 












1 i ! 










0 


Ui T 


O' 










i 








1 




_ 


















; j 


i 








o > j 


t \j 






l ? 1 M 


-r 1 

n 


M ! i 
























i 












I 








*=%.' ; "M 

«cT-» ._ 




s. i 


i 














• -H 






1 




















| 


r- - 

I 


! 




rfrM-l 




1 i I 




! i j 




i 








— \ — 




" 








MM 






. V 
















1 








■ 


















] 








- 


j 


^ u 


,*= ' 






I ! i'l 






_ 




























— 1 






| 




! 









§ 


5 


3 

> 


w u 






ti 


r» 


































— ! 




1 








^ 1 




* 


O Qd 


§ 






> 


i 


















































o 




D? 


































































4 




•-4 






1 




















































~ I 






! 




f 






















| 


i H 






J 






























































1 






















1 


1 








7 










“T 


































~ J 
















1 
































! 


— 






































































































































n 






















' 










i 
















1 






1 














































































1 


















































































; 
















































































































i 




























































































■ 


l j 








l 




i 


~ 


r\ 
















I 




























: 








i 






















■ 


i j 


! 




















j j 


J 































M 










" I | 

| 






1 


1 












I ; 








M_ 


I 

r 




















n 


1 








LI 




1 










(- 




jZflj 




L 


y. 


“T 




M 


tZ 




tj; 


1 














(, 




r 




1 - 




1 M 


_n 


r— 




n- 




I it T i 




M 








I It 






7T T 


If 1 






m 


~ 1 


i i f 


n^r 


r 


i n f 




1 1 1 1 






























































i 










1 










i 
















































i 
































1 
















r 














































































































































































i 
















i 








































































IT 
















1 






















































1 
























| 
















| 


L 




































































































1 






















































































I 




L 




- 


1 














































r 


























j 








i 






t 










~t~ 




























































1 
























































c 












* 






















t" 
















] V 


1 








] 






! 


























J . 




I 








T 
























































7^1 






( 
























i 




i 










1 ^ 










. 1 








L 




L 


M 




1 


r 


■ 




i 


rT 










1 




















L 




[ r r 




i 




i i i 




| 












4) 

Mn 




















1 












r 


r . 

i 






I 


pi 




i i 




_ 


L U MH 


1 


h 3 








1 




















1 1 












M i 


i 




L 
















1 








i 


r 


... j- 


1 1 






i -IT EL4<i_kJ 


- i ! 
; j 


1 

i 


| 








1 


I ; 1 1 






1 










j | | | 


i | T i 






: i j 




h 


r 


lTl-n 




r 


Mi 


r 




rr 






Mil 


1 IT i 


!' i i r 






i 







APPENDIX B LiTEIlATiJilE CITATIONS 



1. Vo a Kannia, "The Impact on Seaplane Floats DstsHIng Landing", 

N.A.C.A., T. N. 321, 1823; 

2. Chu, Wen-Hwa, and Abramson, H. Hormaa, 'Hydrodynamic 
Theories of Ship Slamming-fleview and Extentioa", Journal of 
Ship Research, Vol, 4, No. 4, March 1361; 

3. Sseebehefy, V. G. , and Todd, M, A., "Ship Slamming in Head 
Seas", David Taylor Modal Basin Report 913, February 1855;, 

4. Todd, M. A., "Slamming Due to Part Pitching Motion*", David 
Taylor Model Basin Report c83, January 1855; 

5. Ssebeholy, V, G., and Lum, S. M, Y. , "Modal Sxperimaats on 
The Slamming of a Liberty Ship in Head Seas", David Taylor 
Model B&cin Report 914, February 1055; 

6. Ochi, K&suo, "Model Experiments on Ship Strength and Slamming 
in Regular Waves"', Transactions, SHAME, 135S; 

7. Pierson, J„ D. , "The Penetration of a Fluid Surface by a Wedge", 
Institute of Aeronautical Sciences, Fairchild Fund Preprint No. 
FF-S, July 1950; 

8. Bledsoe, M. D. , "Series Investigation of Slamming Pressures*', 
David Taylor Model Basin Report 1043, December 1958; 

9. Lewis, Frank M. , "The Inertia of the Water Surrounding a 
Vibrating Ship”, Transactions, SHAME, 1929; 

10. Weiabluin, G., and St. Denis, M. , "On The Motion* ii£ Ship# at 

Sea", Transaction*, SHAME, 1950; 

11. Saebehely, V, G,, and Lum, S. M. Y„ , "Model Experlaaents 

on Slamming of a Liberty Ship in Head Seas**, David Taylor 
Model Report 914, February IS 55; 

12. Lamb, Horace, "Hydrodynamics”, Gth Ed. , Dover Publications, 
New York, 1948; 



52 - 



jl JirATl D a mi i „ 'J TL1 a ' - 1 l.v 

, t ■ ‘vcjkl tyihnfO moW • no T :•» *T" .tn**! r._ 7 

,X|S .J5 .T ..a.^.A.W 

r».?C'^n\ > t»o'2fcv>x ‘ ,fl«aa i To>i .£» ,n - •*s3-‘ fact* ,ariO 

&. iMiiKolr s Jtjjft: x3 itelft • it .n_ 1 %i£i 1© *•!***« rTT 

, i 8 % 1 J . 'll > Jw * 4irray6t.fi qi" ' 

**4ft Ki sr !.&«»£* q?r? »,A , JEjfool t .J ,V - 
;!$'l ■**«*»• I .Cl* noqol- nl*& i i»feo* I f. * ,'* 1 fci sC 

b’/uJ ^ saottoM *^2x10*14 tu 1 ! ot *«C v ^ .A ,£i ,fc*L«T 

•8d»x /, v ’ * . f t « *,>".£ M*r.A| %oS^»T 

m zfm, hskj-lL ».7 % .,j . 7 

*.:•■»•#■•.; .' *. jf>a bass:! • q.d* -.J-r- ••' - j *” ixl.-traut?-' oifT 

;SS61 vjt *&< j »m niojetr laboM 

ba« bl*4*''i'T?ls qixi-3 so •Jijoarm* -;c 4 . IrboM ' ,ouuH ,i/b3 
£ * 4 vi t AMAKi t eaois > snail* * t»v? '* xU^-sMt al 

, *■**>«,/ . , : .j • feb-l** « to noi*ant>n 4 *C*~' , 4 ,t ..a©****! 4 ! 

j.rx hi' • !•* * bi-'i . blish%te'Z m ^or»s»it.i U' I© o< ;.Asj4 

-j.’^ vL,-: ,1-^ 

t a**?.*... •<*«* >..'&* ->.UH V ,.Q ,u. 

;?a »tX rc*'J.Tis»'»3C,' ,21*6 1 rs-'<f» ■ .1* *£. lob© W *ob£»T bir- >3 

# */• i ■';<• '7«J , « , W r &tii to m ftxTZ ,,u' abura*! , ji< «d 

5c > ,.tj v !«, cx*r '.il »•_?“ t ,i# ,Jc C>*j ,,%• »r:i.'! iii- iV 

;08SX ,t J/AVli. ,aMoit«u*ur*a*iT *'a»w 

im* ar >* %«u cJT 1 O* t 4 #Y * J 1 * 4 Ji*.I 4 t '•*' * (* J[^*T * -* , 1 ‘ 1 j 

’ *iy» ■ fc , »«#2 b«*H al a.m x&*dlh m lo sxsiuimjsi? ao 

;..c*Si ^s«trid* » >i* bsoq* -< lab ^i.v 

* liicw i <xsvr* , .t-IS iiJ3 83H *n^boTbvb’ ,d<uij 

;i-8i „noY • .1 



13 . 



Pierson, John D. , s ’Cn the Pressure Diatribuxion lor a Wedge 
Penetrating a Fluid Surface", EXT Report No. 3$6 5 Sherman 
M. Fairchild Publication Fund Paper Ho, 19?, Institute of 
Aeronautical Science, New York; 

14. Bisplinghoff, E. 3U , and Doherty, C. S. , "A Two-Dimensional 

Study of th© Impact of Wedges on a Water Surface", Aeroelastic 
and Structures Research, MIT Contract No. adl-SOUl, March 
1050; 

15. Hildebrand, F. B, , "Introduction to Numerical Analysis”, 
McGraw-Hill Book Co. , 1958; 

18. Sokolnikoff, 1. 3., and Redheffer, £. M., ’'Mathematics of 
Physics and Modern Engineering”, McGraw-Hill Book Company, 
1958; 

1?. Borg, S. F., "Some Contributions to the Wedge Water Entry 
Problem", Journal of Engineering Mechanics Division, Pro- 
ceedings ASCJB, Vol, 83, 3MI, April 195?; 

IS. Ssebehely, V. G., "Hydrodynamics of Slamming of Ships”, 

David Taylor Model Basin Report 82S, July 1952; 

19. M i In ® - Thoms on, L, M,, "Theoretical Hydrodynamics 15 , 4th Ed. , 

Macmillan Company, New York, 1980, 



- 55 - 



u -ul r. iu- • '.i« - j is j , — rrTT .SI 

• *-< 1 * ‘'T% »' »1 j i . • t i** >(*•.'* » 

'' " r i 7l !ci - .K 

i± x l . . * .-vi *:*/**♦ A 

1 '* ~ • / 0 n<r*i •*.-» 4 .j .« ... .M 

****** • ?■ ViUft • c-i »- 1W !.• ^ , li UL- tc* vSh./’- 

« A .y , 1 £ W •' f W •..•• •. I* > < .. 



a* J.i. < l rf •*?* ? • •. s 

; - ''i , .»'■** -< 1 

j> A -> * *r-. ... ,.fc *. . ,? i a* , a i 4 #i 

«V*"’< r * ” <,j ‘"l-niiiV,' 1 ir,"vj vi; Tt *»»■>• If ' ' . J* ^ 

t 1 ' 1 , >J,/ 4,: -J.v }. :i *«*'.' .1 ,4 Jica .?! 

■**'* 1 »#'34«. :i *--< -cr - , , n.v.i**} 

i< ?'* . ' ' .If*' * .vi£**a* 

« '•!> ' *'-• ._■/"••*■- l 41 - c' ,.fr .V ;ti 

‘■fe'i . '.i/5* ,<X' K . 4 *. t t *lv. < ia|» i i’ u:ikt 

•-■- . ’ V. r >. T' /' .’ - , , “i ,. \f v{ 

.».'4 J ,> nil <-»k #’< 






•rtf 



thesS924 

Theoretical investigation of slamming lo 




3 2768 002 06017 0 
DUDLEY KNOX LIBRARY 



